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Abstract 

Equilibrium properties of electrons in double-heterojunction AlGaAs/GaAs- 



o 
o 

CO 

o 

I/"") ' /AlGaAs structures are investigated theoretically, using a full self-consistent 

ON ' 



numerical method. The transition from single to bilayer electron systems is 
discussed for the case of double-junction structures with varying distance be- 



tween interfaces. In strong in-plane magnetic fields we observed the same 

o ' 

type of transition accompanied by large effects on the energy spectra. A new 
fenomenon in bilayer systems — separability of the energy dispersion curve 
in parts corresponding to electrons either in first or second layer is analyzed 
in detail. The magnetic field induced variation in the effective distance be- 
tween electron layers is compared to the constant layers separation in standard 
double-well structures. Due to relatively small and soft electrostatic barrier 
between interfaces we expect the van Hove singularity in the 2D density of 
states to be easier detectable than in double- well systems. The cyclotron 
mass, proportional to the density of states and characterizing the electron 
motion in tilted magnetic fields, is calculated as a function of in-plane mag- 
netic field. The separability of the energy spectra resulting in tilted cyclotron 
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orbits is demonstrated on a classical picture of an electron moving in crossed 

electric and magnetic fields. 
PACS numbers: 72.40, 72.20 
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I. INTRODUCTION 



The electronic structure of bilayer 2D electron systems in double quantum wells is in the 
focus of interest since their preparation and intensive theoretical and experimental inves- 
tigation of their properties is now in progress. Main directions of recent research of these 
structures subjected to magnetic fields were reviewed by Eisenstein M. Firstly, in a strong 
perpendicular magnetic field, the double layer systems exhibit new quantum Hall states 
not present in a standard 2D gas because of an extra degree of freedom, the layer index. 
Secondly, in the in-plane magnetic fields, the tunneling between two electron layers can be 
modified to such extent that, at certain critical fields, the interlayer tunneling is completely 
suppressed and electrons move exclusively in one of the wells. The field-induced decoupling 
of electron layers is accompanied by a number of new phenomena which can be studied by 
measuring both the parallel and perpendicular transport and the Shubnikov-de Haas effect. 

A typical double well structure consists of a pair of narrow wells (~ 20 nm) separated by 
an undoped AlGaAs barrier of a similar thickness. Two Si doping layers are usually located 
both below and above the wells and the Schottky gates are deposited on the top/back of 
the sample to control better the concentration of electrons in individual wells. Due to the 
tunneling through the barrier the degeneration of the energy spectrum in the layer index is 
removed, the corresponding single levels split and the lowest bound states of individual wells 
form a symmetric/antisymmetric pair. A small tunneling probability through the hard-wall 
AlGaAs barrier implies the small level splitting and the weak coupling of electron layers. 

Traditionally, the gross features of the electronic structure of a standard double well are 
explained in terms of a simple model which assumes two identical very narrow wells at the 
distance d and a weak overlap of the localized states corresponding to first subbands of both 
wells. The interwell tunneling is described in the tight binding approximation by a single 
parameter determining a small bonding- antibonding energy gap. Then the band spectrum 
is given by the elementary nearly free electron model as presented in classical solid state 
textbooks (see e.g. 0). The model predicts the anticrossing of two circular Fermi lines 
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accompanied by the logarithmic van Hove singularity in the density of states which can be 
interpreted, in the 2D case, in terms of singular effective cyclotron mass. It was utilized to 
study a role of minigaps in Si by Ando || and Mathesson and Higgins Q and successfully 
applied for description of the interwell tunneling (see and references therein) and of the 
electronic structure of double wells in tilted magnetic fields |5|], || and J7|. 

In our theoretical study we shall consider another potentially bilayer systems, the double- 
hetero junction quantum wells formed only by two selectively doped AlGaAs/GaAs interfaces 
with no AlGaAs layer grown between them. Note that in this case an undoped structure 
would result in an exactly rectangular quantum well. In the doped structures, the trans- 
fer of electrons from Si donors into the GaAs leads to the band bending in the GaAs by 
electrostatic forces and, consequently, to the formation of a soft build-in barrier of an ap- 
proximately parabolic shape inside the well. It depends on the well width and the amount 
of the charge wether the built-in barrier is strong enough to separate the quantum well into 
two independent /weakly-coupled parts or if the coupling is so strong that the system can 
be considered as a transient form between a single well and a double well. For simplicity, 
only the symmetrically doped wells will be considered. 

In charged wells, the shape of the barrier is given by the charge distribution. There- 
fore the electronic structure must be determined by the self-consistent solution of coupled 
Schrodinger and Poisson equations. There is a number of such calculations for the case of 
zero magnetic field, see e.g. || and references therein. Effects of the in-plane magnetic field 
on both conduction and valence bands, resulting in a diamagnetic shift of the electron-hole 
recombination energy, were discussed by Oliviera et al. ||. Here the full self-consistent 
calculation will be employed to obtain the electron energy spectra and information about 
the charge redistribution due to the magnetic field. The electronic structure of double- 
heterojunction wells will be studied mainly from the point of view of the possible magnetic 
field induced splitting of a single layer electron system to a bilayer one. 
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II. BAND STRUCTURE CALCULATION 



The standard semi-empirical model working quantitatively for the lowest conduction 
states of GaAs/AlGaAs heterostructures is used to solve the Schrodinger equation in the 
envelope function approximation. The envelope function is assumed to be built from host 
quantum states belonging to a single parabolic band. The effect of the effective mass mis- 
match is completely neglected and the envelope functions of GaAs and AlGaAs are smoothly 
matched at the interface. Due to the translational invariance in the layer plane the momen- 
tum operators p x and p y become good quantum numbers, p x — > hk x and p y — > hk y . If we 
choose a vector potential A in the form A = (B y z, 0, 0), the eigenvalue problem reduces to 
solving the one-dimensional Schrodinger equation with the Hamiltonian 

H = ±pl + ^f(z-z f + V conf (z). (1) 



For more detailed derivation, consistent with our notation, see e.g. [[T(J and references 



therein. The centre Zq of the the magnetic part of the effective potential energy is re- 
lated to the wave vector component k x by zq = hk x /muj, u denotes the cyclotron frequency 
e\B y /m. The confining potential energy of a well with the width d 
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V conf {z) = V b {z) + V sc {z) (2) 

is a sum of the step functions V\,(z) = Vb9(—z — d/2) + V\,Q{z — d/2) corresponding to the 
conduction band discontinuities between AlGaAs and GaAs and of a term V sc (z) describing 
the interaction of an electron with ions and the electron-electron interaction. This term 
should be calculated self-consistently and can be written as 

V sc (z) = V H (z) + V xc (z). (3) 

The Hartree term Vh is the solution to the Poisson equation 

d 2 Vu = \e\g(z) ( 
dz 2 e U 
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and we use an expression calculated by Ruden and Dohler |TT[ in a density-functional for- 
malism for the exchange-correlation term V xc , 

The conduction band offset VJ, and the dielectric constant e enter our calculations as input 
parameters. 

For modulation doped GaAs/AlGaAs heterostructures the total charge density g(z) in 
the equation (§) can be split into parts corresponding to concentrations of electrons, N e (z), 
their parent donors in AlGaAs, N~[(z), and ionized residual acceptors in GaAs, N~(z), 

g (z)=e{N e (z)-N+(z) + N-(z)}. (6) 

We accept a usual approximation of constant impurity concentrations and assume donors 
and acceptors to be ionized within certain finite intervals: N~[(z) = for d/2 + w < \z\ < 
d/2 + w +ld, i-e. symmetrically with respect to the well of the width d, and N~(z) = N a for 
\z\ < d/2, which means that all acceptors inside the well are ionized. Here w denotes the 
spacer thickness and Id is the depletion length of donors which is determined in the course 
of the self-consistent procedure. 



III. RESULTS AND DISCUSSION 

In our calculation we consider four wells with widths d = 25, 40, 50 and 60 nm, close 
to the dimensions of a typical double quantum well described above. The concentrations of 
localized charges are = 2 x 10 18 cm~ 3 , N a = 10 14 cm -3 , the band offset VJ, = 225 meV and 
the dielectric constant e = 12.9. The spacer thickness w = 60 nm yields electron systems 
with the concentration of electrons N e w 3.2 x 10 11 cm -2 . The self-consistent procedure 
was performed in two steps. Firstly, the electronic structure of wells was calculated in zero 
magnetic field, in the second step the magnetic field dependence of subbands was obtained 
starting from intermediate zero-field results. 
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A. Zero magnetic field case 



It follows from the form of Hamiltonian ([!]) that, in the case of zero magnetic field, an 
electron motion in the z direction is completely decoupled from the in-plane components. 
The energy spectrum consists of two quasi-continuous free electron branches, the quadratic 
functions of the wave vector components k x ,k y , and the energy levels corresponding to the 
bound states of a well. Positions of levels on energy scale are presented in figure ffl, together 
with the band diagrams resulting from the first step of the self-consistent procedure. The 
form of the barriers is approximately parabolic for all four wells, the minor deviations from 
this shape may be observed for wider wells. The heights of barriers are very small when 
compared with hard walls of wells formed by the GaAs/AlGaAs interfaces. This indicates 
a stronger coupling between wells than in a standard double-well structure. 

All the systems have one occupied pair of symmetric/antisymmetric subbands, except 
of the 25 nm wide well, where only the subband belonging to the lowest (bonding) state is 
occupied. Note that the origin of the energy scale in figure |l] is located at the Fermi energy 
Ep of the electron system. The energy spectrum of the 60 nm system corresponds to almost 
independent quantum wells, both occupied levels lie below the top of the barrier and the 
splitting of energy levels due to the quantum tunneling is small. The coupling of 2D electron 
layers increases smoothly with the decreasing well width and the energy spectrum of the 
25 nm well reminds rather of a single well than a double well since the energy difference 
between the bonding and the antibonding states is comparable with the difference between 
the energy of the first level and a bottom of the well. 

The electron charge density N e (z) corresponds to well known behaviour of two coupled 
wells as described by simple tight-binding model: there are two maxima of the density 
positioned symmetrically in each well, with the small but finite minimum at z = for 
bonding states and N e (0) = for antibonding states. Of course, the maxima are closer and 
bonding state charge densities higher in narrower double wells. 
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B. Magnetic field dependent band structure 



As seen in figure [I], it is possible to approximate the form of barriers by an expression 
—mQ 2 z 2 /2, where Q is the fitting parameter. Assuming further that the barriers are not 
substantially changed by redistribution of electrons in the magnetic field, the potential en- 
ergy in ([T]) is given by a combination of the fixed electric barrier and the attractive magnetic 
well with the curvature given by B y . The position of minimum, z^, which need not lie inside 
the interval —d/2, +d/2, is determined by k x . For small magnetic fields the potential energy 
inside the well is still a barrier with a smaller curvature m(uj 2 — Q 2 ) and a maximum at 
z = zquo 2 /{uo 2 — Q 2 ). At a critical magnetic field, B y = mQ/\e\, the bottom of the well 
becomes a linear function of z, with the slope dependent on k x , and at higher fields the sign 
of curvature of the potential energy line is reversed and a minimum appears at z instead of 
a maximum. The critical fields for the considered wells B eo = 4.7 T, B 50 = 5.6 T, B m = 6.7 
T and B 2 ^ = 9 T are fully in the range of experimentally available values. This distinguishes 
the double-heterojunction structures from the standard double wells with steep AlGaAs bar- 
riers in which two deep minima are preserved for arbitrary magnetic fields. Consequently, 
the simple tight binding picture of eigenfunctions as symmetric/antisymmetric combinations 
of the eigenstates of individual wells is broken. Above a critical field the lowest eigenstates 
exhibit always a single maximum and thus electrons are more localized by the magnetic 
field around its centroid then in the case of standard structures. Examples of eigenfunctions 
together with 'electro-magnetic' potential energy lines defined by equation flU) are presented 
in figure |2|. The results are plotted for selected values of B y and, to demonstrate their k x 
dependence, also for a minimum and a maximum \k x \ corresponding to occupied states. 

The transition from a single layer to bilayer electron system, induced by the in-plane 
magnetic field, is clearly visible in figure [| Moreover, the effective inter-layer separation 
increases with the magnetic field reaching a value equivalent to the width of the central GaAs 
region, in the strong magnetic field limit. These effects are absent in standard double well 
systems where the distance between 2D electron layers is a constant given by the separation 
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between wells. 

It follows from the above discussion that the in-plane magnetic field does not influence 
the harmonic dependence of subbands on k y at all, but the dependence on k x is changed 
dramatically. Firstly, the energy separation of bottoms of the bonding/antibonding pair 
of subbands increases with growing B y . This process is accompanied by the continuous 
transfer of electrons from the higher occupied subband to the lower one until the second 
subband is completely depleted. Secondly, the curvature of the originally parabolic subbands 
is modified. The antibonding subband remains almost parabolic, similarly as in zero field, 
only its curvature is slightly higher for larger B y . The bonding subband exhibits more varied 
field dependence. The curvature of E(k x ) decreases for k x around k x = 0, at certain value of 
k x reaches zero and then becomes negative. This means that, instead of a minimum which 
corresponds to a zero field, the local maximum develops at k x = 0, accompanied by two new 
minima positioned symmetrically around it. For large enough magnetic fields the energy of 
the maximum becomes greater than the Fermi energy and the spectrum of occupied energies 
splits into two separated parts. This type of behaviour is general for energy spectra of all 
four well widths. In the narrowest well the splitting does not occur in the investigated range 
of fields. Examples of E(k x ) curves for the well of a medium width, d = 40 nm, and several 
selected values of B y are presented in figure |[ 

The electrons with energies close to the Fermi energy Ep are of particular interest since 
they play essential role in the electronic transport properties of the system. The occupied 
and empty states are separated by the Fermi lines in k space. In a zero magnetic field, 
the Fermi contours corresponding to symmetric and antisymmetric subbands have a form 
of two concentric circles. The larger circle describes the bonding and the smaller circle the 
antibonding states. The modification of the shape of E(k x ) curves by the in-plane magnetic 
fields is reflected also in variation of forms of Fermi contours. They acquire a 'peanut' and 
a 'lens' shapes when the sample is subjected to the intermediate magnetic fields. In the 
strong field limit the 'lens' is emptied and the 'peanut' is split into two parts belonging to 
individual wells. This behaviour is illustrated in figure || for the same electron system and 
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the same set of magnetic field values as in figure |^. The deviations from two crossing circles, 
corresponding to a pair of narrow independent wells, are so large that the shape of Fermi 
lines can hardly be described by the above mentioned simple model. 



C. Cyclotron effective mass and real space trajectories 

The separation of electron layers by the in-plane magnetic field may be clearly demon- 
strated if we consider the electronic structure of a system subjected to a tilted magnetic 
field with a strong component B y , parallel to the GaAs/AlGaAs interfaces, and a weak 
component B z oriented perpendicularly. To describe the system, the self-consistent quan- 
tum mechanical calculation of the electron subbands in the presence of B y , presented above, 
is combined with a quasi-classical description of the in plane electron motion under the 



influence of B z , as developed originally by Onsager |12j and Lifshitz JT3[ . 

The semiclassical theory predicts that an electron is driven by the Lorentz force due to 
B z around the Fermi contour with the cyclotron frequency uj c = \e\B z /m c . The cyclotron 
effective mass m c , defined by the above equation, is an important characteristics of each 
Fermi line as a whole and should not be confused with the electron effective mass which is in 
our case a tensor and an anisotropic function of the position on the Fermi line. The explicit 
expression relating m c to the shape of the Fermi contour is 

h 2 r dk 



where dk denotes an element of a length of a Fermi line. In the case of two-dimensional 
electron systems a simple relation holds between the cyclotron effective mass and the density 
of states g corresponding to a single subband: 

irn 

Note that this is the cyclotron effective mass, and not the effective mass, which determines 
e.g. the temperature damping of Shubnikov-de Haas oscillations. 
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Both equations ([?D and ([8]) can be found, at least in an implicit form, in [12| and |13 



already. Since then, they were utilized many times, for application to 2D systems see e.g. [|J, 
but overlooked in the recent publication by Harff et al. [TJj]. Their detailed derivation and 
discussion for the case of combined parallel and perpendicular magnetic fields is presented 



m 
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The magnetic-field dependent cyclotron effective masses resulting from our band struc- 
ture calculation are shown in figure || for all four wells. The logarithmic singularities of 
the cyclotron mass (density of states) curves corresponding to the lowest subbands occur 
at higher fields for narrower wells than for broader ones. The position of a singularity is 
identical with the field value at which a 'peanut' becomes a pair of 'kissing' Fermi contours, 
i.e. when it is to split into two independent parts. Note that this is always above the fields 
Bqq = 4.7 T, B 5Q = 5.6 T and B m = 6.7 T discussed above. The asymmetric shape of van 
Hove singularities somewhat deviates from — ln\EF(B y ) — E \ predicted by the analysis based 
on the quadratic approximation of the energy spectrum near the saddle point. The reason 
for the slow increase of m c on the left side from the singularity, and the rapid decrease on the 
right side , may be due to the large influence of higher order terms in the power expansion of 
E(k x ). For the Fermi contours of the second subband the cyclotron mass always decreases 
for the given concentration of electrons. 

In the real space, electrons move along trajectories having the shape of Fermi contours 
but rotated by an angle ir/2 and multiplied by h/\e\B z . A relation exists between the centre 
of mass (z) of an electron state in a well and the energy spectrum E(k x ) modified by the 
in-plane magnetic field, 

1 ' mu %lo dk x 

Since (z) is a function of k x , i.e. of the position of an electron on the Fermi contour, the 
real space trajectory is slightly tilted with respect to the x — y plane. Travelling around a 
'peanut' Fermi contour with the frequency u c , the electron is transmitted from one well to 
the other when the wavevector k x changes its sign. At k x = the electron appears on the 
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top of the barrier. The separation of a 'peanut' into two parts means that two independent 
real space trajectories for electrons are formed, localized completely in opposite wells. In 
that case, the electron moving around a Fermi contour is always reflected by the barrier and 
remains confined in one of the wells. Hence, the electron system can be considered as divided 
into two disconnected layers. Self-consistently calculated examples of (z) as a function of k x 
are presented in figure ||. 

The quasi-classical model of an electron motion in a double well can help to illustrate 
the origin of localization of electrons in individual wells due to the magnetic field. Figure 
|7| shows real space trajectories of an electron in a well with a parabolic barrier calculated 
classically for the case of two tilted magnetic fields with the same perpendicular component. 
For the field below the singularity, the trajectory overcomes the barrier and an electron 
spends the same time in both wells. Note that, in accord with this statement, the left part 
of the curve (a) lies above and the right part below the projection of the trajectory to the 
x — y plane, (z — 0). The curve (b) is the electron trajectory calculated for the field above 
the singularity (the curves (a) and (b) are offset for clarity). In that case, an electron is kept 
all the time in one well. Just at the magnetic field, when the singularity occurs, electron 
stops at the top of the barrier and stays there infinitely long. It means, in terms utilized 
above, that its cyclotron frequency is zero and the cyclotron mass goes to infinity. 

IV. CONCLUSION 

We have performed self-consistent calculations of the electronic structure of symmetri- 
cally charged double-hetero junction wells in parallel magnetic fields. Their properties are to 
the large extend similar to standard double wells but, on the other hand, there are several 
differences between them which make the double-heterojunctions attractive. 

In both types of structures an in-plane magnetic field substantially modifies the k x depen- 
dence of energy subbands. The following properties result from these field induced changes, 
(i) Electrons are transferred from the antibonding to the bonding band when magnetic field 
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increases, (ii) The Fermi contour topology changes, the originally circular lines acquire a 
'peanut' and a 'lens' shapes and, at certain critical fields, the 'lens' is emptied (first critical 
field) and the 'peanut' splits into two independent lines (second critical field) . (iii) The split- 
ting of the 'peanut' Fermi contour is accompanied by the separation of an electron layer into 
two independent systems, (iv) At the second critical field the density of electron states and 
the cyclotron effective mass exhibit a logarithmic singularity when considered as functions 
of the in-plane field. 

The hard wall barrier formed by the AlGaAs layer means that the coupling between 
electrons in a standard double well is always weak. Therefore, a simple tight binding method 
for description of an electron motion in the z direction, resulting in a nearly free electron 
model for its motion in the x — y plane, can be successfully applied to describe the electronic 
structure. 

This is not true for double- hetero junction systems. Due to the charge present inside 
wells, a soft barrier of approximately parabolic shape is built in each well, with the top in 
its centre. The height of a barrier is very small and the strength of coupling can be tuned 
in large range by proper choice of the concentration of carriers and the width of the well. 
Under these circumstances, the simple model fails to describe the electronic spectra even 
in a qualitative way. Note, that according to our calculations, e.g. the gap between two 
occupied subbans increases under the influence of the field, while the simple model yields a 
constant gap. 

The strong coupling between wells means that stronger magnetic fields are necessary to 
modify the electronic structure of a double-heterojunction. The difference between the first 
and second critical fields can amount several Tesla and, consequently, also the changes of 
the cyclotron effective mass are more robust and probably easier to detect experimentally. 

Nevertheless, the above described differences between two types of structures are mainly 
quantitative. The qualitatively new features of double-heterojunction structures are the 
single layer/bilayer transition and the possible control of the inter-layer separation by in- 
plane magnetic field. These effects might be useful in studies of many-body effects in bilayer 
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electron systems. 

For these properties the double-heterojunction structures represent a very attractive 
alternative to standard double wells. 
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FIGURES 



FIG. 1. Band diagrams and energy levels of double- heteroj unction wells at a zero magnetic field. 

FIG. 2. Wave functions and 'electro-magnetic' band diagrams of a double- heteroj unction well 
subject to in-plane magnetic fields. Pairs of curves characterize the electrons with a minimum 
(dashed) and a maximum (full) occupied \k x \. The width of the well is 40 nm. 

FIG. 3. Occupied energy subbands of a double- heteroj unction well, d = 40 nm, subject to 
in-plane magnetic fields as functions of k x . 

FIG. 4. Fermi contours of a double- heteroj unction well, d = 40 nm, subject to in-plane magnetic 
fields. 

FIG. 5. Cyclotron effective masses of four wells as functions of an in-plane magnetic field. 

FIG. 6. Centroids of electron states of a double- heteroj unction well, d = 40 nm, subject to 
in-plane magnetic fields. 

FIG. 7. Classically calculated trajectories of an electron in a double well subject to tilted 
magnetic field. The upper curve (a) corresponds to B y below the singularity and the lower curve 
(b) to By above this value. 



16 



